Measuring the size of a Schrodinger cat state 
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We propose a measure for the "size" of a Schrodinger cat state, i.e. a quantum superposition of 
two many-body states with (supposedly) macroscopically distinct properties, by counting how many 
single-particle operations are needed to map one state onto the other. This definition gives sensible 
results for simple, analytically tractable cases and is consistent with a previous definition restricted 
to Greenberger-Horne-Zeilinger-like states. We apply our measure to the experimentally relevant, 
nontrivial example of a superconducting three-junction flux qubit put into a superposition of left- 
and right-circulating supercurrent states and find this Schrodinger cat to be surprisingly small. 



Introduction. - In his landmark 1935 paper 
Schrodinger introduced the notion of entanglement, 
and immediately pointed out its implications for 
measurement-like setups, where a microscopic quantum 
superposition may be transferred into a superposition of 
two "macroscopically distinct" states. His metaphor of a 
cat being in a superposition of "dead" and "alive" , ini- 
tially designed just to reveal the bizarre nature of quan- 
tum mechanics, nowadays serves as a namesake and in- 
spiration for a whole generation of experiments designed 
to test the potential limits of quantum mechanics in the 
direction of the transition to the "macroscopic" world, 
as well as to display the experimentalists' prowess in de- 
veloping applications requiring the production of fragile 
superpositions involving many particles. Recent experi- 
ments or proposals of this kind include systems as diverse 
as Rydberg atoms in microwave cavities (2], supercon- 
ducting circuits [111, H, H I3J3 , optomechanical [i, 
and nanoelectromechanical systems, molecule inter- 
ferometers fl2|. magnetic biomolecules , and atom op- 
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The obvious question of just how many particles are in- 
volved in such a superposition has not found any general 
answer so far |0| , and discussions of this point in relation 
to existing experiments have often remained qualitative. 
While the number of atoms participating in a macro- 
scopic superposition of a Ceo molecule being at either 
one of two positions separated by more than its diameter 
is obviously sixty, the mere presence of a large number of 
particles in the system is not sufficient in itself. This is 
demonstrated clearly by the example of a single electron 
being shared by two atoms in a dimer, atop the back- 
ground of a large number of "spectator electrons" in the 
atoms' core shells. Therefore, obtaining a general mea- 
sure for the "size" of a superposition of two many-body 
states is nontrivial, especially for systems such as su- 
perconducting circuits, where the relevant superimposed 
many-body states are not spatially separated. 

Before proposing our solution to this challenge, we 
state right away that certainly more than one reason- 
able definition is conceivable, depending on which fea- 
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Figure 1: (Color online) (a) Example of normal-state per- 
sistent currents mentioned in the text, where D = 3 single- 
particle operations are needed to turn state \A} into \B). (b) 
Hilbert spaces Hd generated by repeated application of single- 
particle operators, (c) Probability distribution Pe{D — d) 
for the distance D between two BEC states or between the 
two components of a generalized GHZ state, as a function of 
the angle between the corresponding single-particle states, for 
N — IQ particles, see Eq. (0J. 



tures of the state are deemed important for the superpo- 
sition to be called "macroscopic". Previous approaches 
can be roughly grouped into two classes: Measures of 
the first kind involve considering some judiciously chosen 
observable, evaluating the difference between its expecta- 
tion values for the two superimposed states, and express- 
ing the result in some appropriate "microscopic units" 
or in units of the spread of the individual wave 
packets Several recent experiments have produced 

Schrodinger cats that, by those measures, are remarkably 
"fat". For example, for the experiments in Delft0| and 
the clock- and counterclockwise circulating su- 
percurrents, whose superposition was studied, were in the 
micro- Ampere range, leading to a difference of lO^/^s or 
even lO^^/xs in the magnetic moments, respectively. 

Measures of the second kind, in contrast, try to infer 
how many particles are effectively involved in those su- 
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perpositions, which will be the focus of our paper. This 
category comprises Leggett's "disconnectivitY;" 
and the measure of Diir, Simon and Cirac [l^ (DSC). 
The latter applies to a class of generalized Greenberger- 
Horne-Zeilinger (GHZ) states, which it compares to stan- 
dard GHZ states in terms of susceptibility to decoherence 
and entanglement content. 

In this paper, we propose a general definition for the 
size of a Schrodinger cat, or the effective distance between 
its two constituent many-body states, that is based on 
asking the following question: "How many single-particle 
operations are needed (on average) to map one of the two 
states to the other?" (see Fig. 1 a) 

We will make this definition precise using the language 
of second quantization and show that in simple analyt- 
ically tractable cases it agrees with reasonable expecta- 
tions and with the measure of DSC IT§| . After analyzing 
the general features of our measure, we evaluate it nu- 
merically for a superconducting three-junction flux qubit, 
whose eigenstates we find by exact numerical diagonal- 
ization. The resulting size turns out to be surprisingly 
small, which we attribute to the fact that repeated ap- 
plications of single-particle operators quickly produce a 
very large Hilbert space, in which the "target" many- 
body state can be represented accurately. 

Definition of the measure. - We start with a simple ex- 
ample. Consider a clean, ballistic, single-channel metallic 
ring capable of supporting normal-state persistent cur- 
rents of electrons. Suppose it is put in a superposition of 
two Slater determinants, |^) and \B), which differ only 
in the number of right- and left-moving electrons (Fig. 
n a) . The number of particles effectively participating in 
this superposition is clearly equal to the number of parti- 
cles that have to be converted from right- to left-movers, 
in order to turn one of these many-body states into the 
other. This is identical to the number D of single-particle 
operations that have to be applied to effect this change: 

I5)«nf=i4'% 1^)- 

When turning this into a general definition, we have to 
realize that the "target" state \B) might be a superposi- 
tion of components that can be created from \A) by ap- 
plying different numbers d of single-particle operations. 
In that case, we will end up with a probability distribu- 
tion P{D = d), defined as the weights of these compo- 
nents, for the "distance" D between \A) and \B) to equal 
d. Furthermore, repeated application of single-particle 
operations may lead to a state that could have been cre- 
ated by a smaller number of such operations (such as 
\B) cx clck'cliCk \A) cx \A) if Uk — 1 and Uk' = 0). This 
has to be taken care of by projecting out the states that 
have been reached already. 

The general definition (whether for fermions or bosons) 
is obtained by starting from the Hilbert space Tio = 
span{|A)}, and applying iteratively the following scheme, 
to generate spaces Tii, 7^2, • • ■ (Fig. fflb): Given a Hilbert 



space Hd, apply all possible single-particle operators c]ci 
to all of its vectors. Consider the span of the resulting 
vectors and subtract the orthogonal projection on all pre- 
vious Hilbert spaces, Ho©Hi© . . .©Hd, thereby generat- 
ing Hd+i- This scheme is guaranteed to produce all vec- 
tors that can be generated from \A) by the time-evolution 
of an arbitrary (interacting, possibly time-dependent, but 
particle-conserving) Hamiltonian. Thus, we can repre- 
sent the "target" state \B) (which is assumed to have 
the same particle number) as a superposition 

oo 

\B)=Y,XM (1) 

d=0 

of orthonormalized vectors \vd) G Ti-d- The amplitudes 
Xd (defined up to a phase) produce the probability dis- 
tribution for the distance D from \A) to \B), 

P{D = d) = \Xd\\ (2) 

from which an average distance D may be obtained. 

Generalized GHZ states. - Before discussing general 
features, let us consider an important example, namely 
a superposition of two non-interacting pure Bose con- 
densates, \A) and \B), with N particles being simulta- 
neously in the single-particle state \a) or respec- 
tively, where {a\P) — cos 6*. We can write the two 
BEG many-body states as \A) = (7V!)-i/2(cl)^ |Vac) 
and \B) = (iV!)-i/2 (cos 6lcl + sm9cl)^ \Yac), with c\ 
creating a particle in state |a), and creating a par- 
ticle in the state that defines the orthogonal direction in 
span{|a) , (we have dropped a potentially present, 
but irrelevant global phase). Expanding \B), we find: 

1^) - 7^ E ( d ) (sin ^^4)' (cosset)'"'' I Vac) . 

(3) 

Then \vd) = [dl{N ~ d)!]-i/2 hlj (cU |Vac) is a 
normalized state that can be reached from \A) in ex- 
actly d applications of the single-particle operator c\ci, 
i.e. \vd) G Ti-d. Thus, \B) may be represented in the form 
(P), with coefficients 




The resulting distribution Pg{D = d) = \Xd\'^ is binomial 
(Fig. He), with probability p = sin^ 6 = 1~ |(a|/3)|2, 
and thus the average distance turns out to be 5 = Np. 
It will be maximal, D ~ N, for two orthogonal single- 
particle states, as expected. This example can be tran- 
scribed into spin-language, by considering the states 
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\A) = It)^ and \B) = (cos 61 |T) +sin 6* In that case, 

we have to adapt our approach by defining single-spin op- 
erators as the single-particle operations, and replace c\ci 



(a) aEj, aC 



(C) 



by X^jLi '5'-^ ■ Straightforward algebra shows the results 
for P{D = d) and D to be the same. Comparing to DSC 
[lij i. where such generalized GHZ-states were analyzed, 
we find that our result agrees precisely with theirs, for 
this special class of states, to which the analysis of DSC 
was restricted. 

General features. - We can prove that the Hilbert 
spaces Tid thus constructed do not depend on the choice 
of the single-particle basis used to define the operators 
Consider an arbitrary unitary change of basis, 
'YlijUijCj. Starting from an arbitrary vector 

we want to show that spanjc^^c^ \v)} = span{c|^Cj 
(where i,j range over the basis). Indeed, any vector \w) 
from the Hilbert space on the left-hand-side can be writ- 
ten as \w) = J2i,^j,^,^jfJ-z',j'U*,^Uj'jclcj\v), which is an 
element of the right-hand-side (and vice versa). As a 
result, no particular basis (e.g. position) is singled out. 

We can prove as well that the distance is symmet- 
ric under interchange of \A) and \B) for an important 
class of states, namely those connected by time-reversal 
(such as left- and right-going current states considered 
below). With respect to a position basis of real- valued 
wave functions, this means \A) — \B)* . In that case, 
since the single-particle operators can be chosen real- 
valued, we have Tif^^ ~ O^d^^)* ^ rnaking the weights 
P^-*^{D = d) and P^^^{D = d) equal. The example 
treated above can also be expressed in this way, by an 
appropriate change of basis, with \A/ B) cx (cos(|)c| ± 

isin(l)cj)^ |Vac). For other, non-symmetric pairs of 
states this is not true any longer, i.e P^~*^ 

can become different from P^^^. An extreme exam- 
ple is provided by the states |^) = ^(|7V, 0) + |0, iV)) 
and \B) — jA^ — 1,1), for A^ bosons on two islands 
(with 1^1,712) denoting the particle numbers). Here 
P^^b[d = 1) = 1, but P^^^{D = 1) < 1, with 
P^^\D = iV - 1) 7^ 0. In the following, we will re- 
strict our discussion to time-reversed pairs of states. 

Application to superconducting circuits. - A super- 
conducting circuit such as a Cooper pair box or a flux 
qubit device can be viewed as a collection of metallic 
islands between which Cooper pairs are allowed to tun- 
nel coherently through Josephson junctions. Adopting 
a bosonic description, we would describe tunneling by 
a term cjcj, where Cj annihilates a Cooper pair on is- 
land j. However, as the total "background" number 
of Cooper pairs n on any island is very large and ulti- 
mately drops out of the calculation, the more convenient 
(and standard) approach is to consider instead opera- 
tors 6"*"^^ = X^ra W ~ ('^Ij that reduce the number of 
Cooper pairs on island j by exactly one. Then, the tun- 
neling term becomes n~^c\cj ^ e^ivi-Vj) ^ while the total 
electrostatic energy may be expressed by the number op- 
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Figure 2: (Color online) (a) Circuit diagram of the three- 
junction flux qubit. (b) Equivalent representation in the 
charge basis, (c) Energy-level diagram for Ej/Ec = 20 and 
Of = 1, as a function of magnetic frustration. At / = 0.5, the 
ground and first excited state are superpositions of left- and 
right-going current states, [±/), the states between which we 
calculate the "distance" D. The current distribution in the 
ground state is displayed in the inset. 



erators fij that count the number of excess Cooper pairs 
on each island. These two types of operators define the 
single-particle operators needed in our approach. 

Let us now apply the measure defined above to a par- 
ticular, experimentally relevant case, the three-junction 
flux qubit that has been developed in Delft 0, H^ . 
Three superconducting islands are connected by tunnel 
junctions (Fig. [21, where the tunneling amplitude is 
given by the Josephson coupling Ej^ and the charging 
energy Ec = e^/2C is determined by the capacitance C 
of the junctions. One of the junctions is smaller by a 
factor of a, introducing an asymmetry that is important 
for the operation of the device as a qubit. The tunneling 
term in the Hamiltonian is given by 



H 



^ . . 

where the externally applied magnetic flux $ = /<I>o is 

expressed in units of the flux quantum $0 = h/(2\e\) to 
define the frustration / that enters the extra tunneling 
phase 6 = 27r/. The charging Hamiltonian is 



(Qi - Qa)^ 
2 + 1/ a 



(6) 



with Qj = 2enj and the restriction X]j=i Qj ~ 0- 
simplicity, we have neglected the small effects of the self- 
inductance and external gate electrodes. 

At / = 0.5, the classical left- and right-going current 
states are degenerate in energy, and quantum tunnel- 
ing (via the charging term) leads to an avoided crossing, 
with the ground and first excited state becoming super- 
positions of the two current states. We diagonalize the 
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Figure 3: (Color online) (a) Average many-body distance D 
between the left- and right-going current states forming the 
ground state of a three-junction flux qubit at / = 0.5, plotted 
as a function of Ej/Ec, for various asymmetry parameters 
a. (b) Corresponding probability distribution for a — 0.8. 
(c) Magnitude / of the average current in the two current 
states, and average charge fluctuation 5N in the ground state 
(symbols as in (a)). 



current operator / = —dH/d^ in the two-dimensional 
subspace of the ground- and first excited states, which 
results in eigenvalues ±J belonging to the two counter- 
propagating current states |±/). Whenever the excited 
and the ground state are well removed from higher lying 
levels (as should be the case in a flux qubit), an equiva- 
lent way of finding |±/) is to write the ground state as a 
superposition of current operator eigenstates in the full 
Hilbert space, and keeping only the contributions with 
positive or negative current eigenvalues, respectively (as 
indicated in the inset of Fig. 12) ■ The distance D between 
the states |±/) then provides a measure of how "macro- 
scopic" the ground (or excited) state superposition is, in 
the sense of the approach outlined above. 

Our calculations have been performed in the charge ba- 
sis, by truncating the Hilbert space to (2An -I- 1)^ states 
|ni, 712, n^) with ni^2 = —An ... -I- An (and — — ni — 
n2). Exact numerical diagonalization of H = Hj + H^h 
yields the ground state and the first excited state, and, 
from them, the current states |±/), as explained above. 
The approach is then implemented by applying itera- 
tively all possible single-particle operators (represented 
as matrices), starting from \ A) = The target state 

\B) = |— /) is represented as a superposition in the 
resulting Hilbert spaces Hd, which yields the weights 
P{D = d). 

FigureOshows the resulting average distance D (calcu- 
lated with An = 6). The fact that D > 1 is a consequence 
of defining the two states \A) and \B) as the eigenstates 
of the hermitean current operator, which makes them or- 
thogonal by default, thus P{D = 0) = 0. Ata = l, 
the monotonic rise of D with Ej/Ec is expected, as a 



larger Ej/Ec allows the charges on each island to fiuc- 
tuate more strongly, implying that more Cooper pairs 
can effectively contribute to the current states. The non- 
monotonic dependence on Ej/Eq for a < 1 was unex- 
pected, but is likely due to the fact that smaller values 
of a tend to make the two counterpropagating current 
states no longer a "good" basis (the ring is broken for 
a — 0). In Fig. 13 (c), we have plotted both the ex- 
pectation value of the current operator in one of the 
two superimposed states, as well as the average parti- 
cle number fiuctuation 5N in the ground state, where 

5N^ = |Ej=i (("j - (%>)^)- Evidently, neither of 
these quantities can be directly correlated to the aver- 
age distance Z?, apart from the general trend for all of 
them to usually increase with increasing Ej/Eq- 

What is initially surprising is the fact that the distance 
remains small, although the examples discussed earlier 
clearly show that much larger distances may be reached 
in principle when applying our measure. In contrast, the 
"disconnectivity" for the Delft setup was estimated [l^ 
to be on the order of 10^, although a rigorous calculation 
seems very hard to do. Two reasons underly our finding 
for the fiux qubit: First, it appears that the fiux qubit 
considered here is really not that far from the Cooper pair 
box. In the Cooper pair box 20], only a single Cooper 
pair tunnels between two superconducting islands, yield- 
ing D — \. In fact, allowing only for a small charge 
fiuctuation (e.g. An ~ 4) on each island of the fiux qubit 
system is sufficient to reproduce the exact low-lying en- 
ergy levels of this Hamiltonian to high accuracy for the 
parameter range considered here, since the charge fluctu- 
ations grow only slowly with Ej/Ec, as observed in Fig. 
131(c) (SN ~ (Ej/Ec)^/'^ at large Ej/Ec). This means 
from the outset that very large values for D may not be 
expected. Second, when analyzing the structure of the 
generated Hilbert spaces Hd, it becomes clear that the 
dimensions of those spaces grow very fast with d, due to 
the large number of combinations of different single par- 
ticle operators that are applied. For that reason, it turns 
out that the "target state" \B) — |— /) can accurately be 
represented as a superposition of vectors lying within the 
first few of those spaces, yielding a rather small distance 
D. 

Outlook. - Our measure for the "size" of Schrodinger 
cat states can be applied, in principle, to any superposi- 
tion of many-body states with constant particle number. 
Future challenges include the extension to states without 
a fixed particle number and the comparison to other mea- 
sures, besides the DSC result 0. In those cases where 
different particles couple to independent environments 
(as was assumed in DSC), our measure is expected to 
be an indication of the decoherence rate with which the 
corresponding superposition is destroyed, and it would 
be interesting to verify this in specific cases. Finally, we 
note that applications to many other physical systems 
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arc in principle straightforward, though the fast growth 
of Hilbert spaces may represent a significant hurdle for 
the direct numerical approach used here, and more effi- 
cient approximations would be desirable. 
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